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Abstract 


Influence of an Integral Splitter Plate and an An- 
nular Fin on Heat and Momentum Transport for 
the Case of Flow Past a Tube 

The longitudinal finning of a circular tube is assumed to be in a configuration 
where the fin is attached at the back of a circular tube The longitudinal fin, 
attached at the back of the tube, is called integral splitter plate The fiow and 
heat transfer results are presented for three different lengths (L=2H, 3H and 4H) 
of the splitter plate and three different Reynolds numbers (500, 1000 and 1500) 
A splitter plate or longitudinal fin connected integrally to a circular tube (behind 
the tube) plays an important role m reducing the size of the wake zone The 
flow physics of such a system becomes quite interesting The role of the splitter 
plate m heat transfer enhancement is an important aspect to be investigated 
The wake of a circular tube without a splitter plate removes heat from the tube 
surface more effectively as compared to the present case of tube with splitter 
plate However the combined effect of conduction and convection m the latter 
case IS certainly greater than that of the former one 

We focus our attention on a section where only one annular fin has been con- 
sidered The geometry of interest is basically a representative section whereby 
the length of the tube is equal to the pitch of the annular fins Appropriate grid 
is generated and computations are carried out for Ren = 1000 The effect of 
annular fin on heat transfer and flow characteristics is to be compared with the 
heat transfer and flow characteristics for the case of a circular tube The results 
are presented for annular fin of outer diameter 1 5 times the diameter of the tube 
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Chapter 1 
Introduction 

1.1 Problem Description 

Fm-tubes are commonly used in many gas-liquid cross-flow heat exchangers In 
such an arrangement, the cold stream (gas) generally flows across outside the 
tubes and the hot fluid (often it is liquid) flows inside the tubes The fins provide 
the extra heat transfer surface area for the cold fluid Even with the extended 
surfaces, the dommant thermal resistance m air cooled condensers is on the air 
side, since condensation (phase-change heat transfer) takes place inside the tubes 
Therefore, in order to achieve signiflcant heat transfer enhancement, strategies 
must be developed, that would result m increased heat transfer coefficients on 
the fin and the tube outer surfaces without a large increase m pressure drop. 

Figure 1 1 shows a schematic diagram of the core region of a fin-tube heat ex- 
changer Numerical investigations m this field have been carried out by Biswas 
et al (1994), where they have tried to achieve the enhancement of heat trans- 
fer from the fin surfaces by inducing longitudinal streamwise vortices in the flow 
field These vortices were generated with the help of delta-winglet type vortex 
generators on the flat surface (Figure 1 2) The longitudinal vortices are developed 
along the side edge of the delta-winglets due to pressure difference between the 
front surface facing the flow and the back surface Sudi vortices are also called 
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streamwise vortices because they have axes parallel to the flow direction Once 
they get developed, they interact with an otherwise a two-dimensional boundary 
layer and produce a three-dimensional swirling flow that mixes near wall fluid 
with the free stream This is a type of kinematic mixing and strongly enhances 
the entrainment of fluid from periphery to the core region of the flow field It 
results m disruption of thermal boundary layer which causes enhancement of heat 
transfer The additional pressure losses happen to be modest because the form 
drag for such winglet-type slender bodies is low 

In the present work numerical simulations have been carried out for flow and 
heat transfer in fin-tube cross-flow heat exchangers We have considered two dif- 
ferent flow configurations In one case longitudinal fin is attached at the back 
of a tube and in another annular fins are attached on the tube The grid-mesh 
used is of body-fitting type employing finite volume algorithm of Eswaran and 
Prakash (1998) 


1.2 The Scope of Present Work 

In order to achieve the objective outlined obove, a numerical investigation has 
been carried out The complete Navier-Stokes equations together with the energy 
equation are solved in a rectangular channel with a built-in circular tube having 
a longitudinal fin, i e , Splitter Plate and built-in circular tube having a circum- 
ferential fin, i e , Annular Fin A detailed analysis of the flow structure together 
with heat transfer characteristics in such modules have been performed. 
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1.3 Outline of the Thesis 

In Chapter-1, the genesis of the problem has been discussed The Chapter-2 of 
the thesis provides a review of the literature relevant to understanding of the ba- 
sic mechanism involved m augmentation of heat transfer, flow and the algorithms 
related to solving the Navier-Stokes equations The mathematical formulation 
of the problem for simulation is presented in Chapter-3 In this Chapter, the 
flow geometry, the governing equations, boundary conditions and grid generation 
techniques are discussed Chapter-4 discusses the solution algorithm and all the 
basic description of the finite volume formulation m detail Chapter-5 is divided 
into two sections In the first section a detailed study of flow physics and asso- 
ciated heat transfer are discussed for the case of a circular tube with an integral 
splitter plate and m the second section the same is discussed for the case of a 
circular tube with an annular fin Chapter-6 includes the concluding remarks and 
the scope for future research 


Liquid 
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Gas 


Figure 1 1 Schematic diagram of core region of a fin-tube heat exchanger 



Winglet pair 


Figure 1 2 Delta-Winglet type vortex generators on flat surface 






Figure 1 3 Three-dimentional Computational Domain for Splitter Plate 



Figure 1 4. Three-dimentional Computational Domain for Annular Fm 



Chapter 2 

Literature Survey 


2.1 Introduction 

Considerable amount of experimental, anal3dical and computational research has 
been earned out on enhancement of heat transfer The enhancement of heat 
transfer through manipulation of surface geometry has been analyzed by many 
researchers and practitioners In order to analyze the flow structure and heat 
transfer m such applications, a detailed computational study is needed The 
flow regime may possess either laminar or turbulent characteristics Numerical 
solutions for laminar flows and heat transfer have been well developed by various 
researchers but the corresponding investigations for turbulent flows involve more 
rigorous considerations This chapter outlines a summary of articles available m 
the areas closely related to the present investigations The review of literature 
is divided in to three different sections In the first section, an overview of work 
done by various researchers in the area of enhancement of heat transfer and flow 
structures are presented There are three subsections m the first section The first 
subsection focuses on studies pertaimng to the Vortex Generators The second 
one describes the work related to integral Splitter Plates and the third one is on 
Annular Fm The second section focuses on the numerical methods with regard 
to the task of computing flow fields m complex geometnes Finally the benefits of 
current problem is established from the available knowledge of existing literature 
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2.2 Enhancement of Heat Transfer and Effect 
on Flow Characteristics 

The primary interest of our present work is to enhance the heat transfer in the gas 
side of fin-tube heat exchangers keeping the pressure-drop penalty at a modest 
level With this intent, we would like to study different investigations related 
to augmentation of heat transfer suitable for the applications enumerated earlier 
Augmentation of heat transfer is of special interest in channel flows where the rate 
of heat transfer between the fluid and channel walls deteriorates as the boundary 
layer grows on the channel walls and the flow tends to become fully developed 
Protrusions can be mounted on these channel walls in order to disrupt the growth 
of boundary layer and thereby enhance the heat transfer between the flowing fluid 
and channel walls. Two relevant applications using such flow configurations are 
the heat transfer between the gas and the fin m the case of gas-liquid fin-tube 
crossflow heat exchangers and the heat transfer between the flowing fluid and 
plates m the case of fin-plate heat exchangers 

2.2.1 Vortex Generators 

The use of a multi-layered surface geometry for plates can disturb evolution to- 
wards a fully developed flow A detailed performance data for louvered fin surfaces 
IS provided by Achaichia and Cowell (1988) In the louvered fins, enhancement 
in heat transfer is obtained at the cost of a high pressure drop In order to cir- 
cumvent this diflS-Culty, protrutions in the form of slender delta-wmglets can be 
mounted on the fin surfaces (Figure 1 2) As shown, the base of the wing remains 
attached to the fin and the apex faces the incoming stream with an angle of 
attack The longitudinal vortices are generated along the side edge of the wing- 
shaped vortex generator due to the pressure difference between the front surface 
facing the flow and the back surface These longitudinal vortices, generated by 
the vortex generators, can be made to disturb the growth of boundary layer in a 
dhannel by exchanging the fluid from the near- wall-region with the channel-core- 
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region and thus they can serve to enhance the heat transfer rate while producing 
less of a pressure drop Use of longitudinal vortices for boundary control is well 
known (Pearcey, 1961) and the vortex generators are used in commercial airplanes 
for this purpose 

Detailed investigation of the enhancement potential of vortex generators have 
been carried out m the last two decades (Fiebig et al (1986), Turk and Junkhan 
(1996), Eibeck and Eaton (1986), Fiebig et al (1989), Amon (1989), Dong 
(1989), Zhang (1989), Tigglebeck (1992), Yanagihara and Tom (1990), Fiebig 
et al (1990), Fiebig et al (1991), Sanchez (1989), Valencia (1992), Biswas et 
al (1994) and Valencia (1996)) The idea of using vortex generators m heat ex- 
changer IS being attempted by the industry Transverse vortex generators, such 
as ribs and corrugations m channel generate vortices whose axes are mainly trans- 
verse to the primary flow direction, whereas, longitudinal vortex generators, such 
as delta wings and rectangular wings generate vortices having axes along the pri- 
mary flow direction The experimental and numerical investigations conducted so 
far show that the longitudinal vortex generators are better options as compared 
to the transverse vortex generators for improvement of heat transfer surfaces in 
the heat exchangers when both heat transfer enhancement and flow loss are taken 
into account 


2.2.2 Splitter Plate 

Experiments reported by a number of authors have shown that the charactenstics 
of the separated wake downstream from a bluff cylinder can be greatly affected 
when a splitter plate is placed along the centerline of the wake. When a splitter 
plate IS placed downstream of a circular cylinder in a crossflow at Reynolds num- 
bers (at which a regular vortex street is shed from the plane cylinder), the vortex 
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shedding may be altered or even suppressed and the drag force experienced by 
the cylinder may be affected 

Roshko (1954) did experiments with splitter plate in the wake of a circular cylin- 
der m a two-dimentional cross flow at a Reynolds number of 1 45 x 10^ He found 
that a splitter plate of length 5D (where D is the diameter of the cylmder)in 
contact with cylinder inhibited the periodic vortex formation and caused the re- 
duction in pressure drag experienced by the cylinder A splitter plate of length 
1 14P m contact with the cylinder did not inhibit the vortex formation, but it 
caused an increase in the base pressure and a reduction in the Strouhal number 
as compared with the cylinder alone Roshko observed a complicated and inter- 
esting sequence of changes in the base pressure on the cylinder surface and in 
Strouhal number when he moved the shorter splitter plate downstream, leaving 
a gap between it and the cylinder 

Gerrard (1966) measured the frequency of vortex sheddmg from a circular cylin- 
der in a cross-flow at Reynolds number of 2 x 10'* with splitter plate of different 
lengths upto a maximum of 2D attached to the cylinder As the splitter plate 
length was increased, the Strouhal number was found to decrease to a minimum 
of St = 0.18 for a plate length of approximately D Further increase in length 
beyond LfD = 1 resulted in a progressive increase in sheddmg frequency 

Apelt and Isaacs (1970) found that short splitter plates attached to a circular 
cylinder in a cross-flow at a Reynolds number of 1 58 * 10* caused large reduction 
in drag With a splitter plate of length 0 125D the pressure drag was reduced 
to 83% of the value measured for the plain cylinder For a sphtter plate length 
of D the pressure drag was reduced to 68% of the plain cylinder value, but only 
relatively small further changes in drag resulted when the splitter plate length 
was increased in steps up to 8D 
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Apelt et al (1973) studied experimentally the effect of splitter plates on the 
flow past a circular cylinder in the range of 10^ < i2e < 5 x 10^ They concluded 
that the splitter plates reduce the drag markedly by stabilizing the separation 
points and produce a wake narrower than that for the case of cylinder alone The 
splitter plates raise the base pressure by as much as 50% and affect the Strouhl 
number to a lesser degree While Cjd and St have minimum values at about 
L/D = 1, the base pressure, Cp was maximum 

Sparrow and Kang (1985) conducted experimental investigations on longitudinally- 
finned cross-flow tube banks and studied heat transfer and pressure drop char- 
acteristics They showed that a high degree of heat transfer enhancement was 
obtained by finning the tubes The finning related enhancements were compared 
with those attainable by increased diameter unfinned tubes Finning was found 
especially advantageous when the comparison was made at a fixed pressure drop 
They considered longitudinal fins at the front, at the rear or both at front and 
rear The rearward fins, which are basically wake splitters, produced a lower 
pressure drop than all other arrangements The experiments were carried out in 
the Reynolds number (based on tube diameter) range of 1200 to 8600 

Mansingh and Oosthuizen(1990) earned out experimental study for the effect 
of splitter plates on the wake behind a two-dimensional bluff body (a rectangu- 
lar cylinder) with fixed separation points at low Reynolds numbers between 350 
and 1150 The results indicated that splitter plates alter the manner of vortex 
formation in the wake causing a decrease in shedding frequency, an mcrease in 
base pressure and a reduction in the overall drag by nearly 50% 

Anderson and Szewczyk (1997) conducted experiments in an atmospheric wmd 
tunnel to study the near wake of a circular cylinder at subcritical Reynolds num- 
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bers between 2700 to 46000 A base mounted splitter plate allowed for the mod- 
ification of the formation region characteristics without disrupting the normal 
Karman shedding Their results provided an explanation for the non-lmearity 
in the relationship between shedding frequency and the splitter plate length In 
addition to the nominal 2-D configurations, a sinuous trailing edge Splitter Plate, 
cylinder taper and shear flow were incorporated to study the three-dimensionality 


2.2.3 Annular Fin 

Transverse annular-finned tubes are widely used as a means of enhancing cov- 
ective heat transfer in a crossflow The flow around an annular-finned circular 
tube in a crossflow exhibits extremely complex three-dimentional flow character- 
istics Hyung et al (1995) carried out experiments for local mass transfer from 
a circular tube, which is situated between two annular fins m a crossflow, at a 
over moderate Reynolds numbers (33 x 10^ < Rcd < 80 x 10®) They employed 
naphthalene sublimation technique to measure the local convective mass transfer 
around the circular tube The ratio of axial gap distance between two axial fins 
(L) to the radial protrusion length of the annular fin (jH’)was an important ex- 
perimental parameter and its range was 0 < L/H <0 40 They found that for a 
small gap distance {L/H = 0.05), naphthalene sublimation is intense in the front 
portion of the circular tube, especially near the separation region. When the gap 
distance is moderate {L/H = 0 15), the overall character of sublimation shows a 
more organized pattern. The influence of horseshoe vortices is substantial For 
a large distance {L/H = 0 40), the impact of horseshoe vortices on the comer 
junction in the vicinity of the endwall is substantial 

Sparrow and Chastain (1986) conducted wind-tunnel experiments to investigate 
the sensitivity of the heat transfer coefficient for an annular fin to small angles 
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of attack (—2° < a < 2°, Q!=angle of attack) They found that the heat transfei 
coefficient m forward zone of the tube was highly sensitive to the angle of attack, 
increasing by about 50% on the face of the fin where leading-edge separation was 
activated as the angle was varied from —1° to 2° The increase in heat transfer 
coefficient at the side of the tube and behind the tube were about 10% Their 
experiments showed that, for any angle of attack in the investigated range, the 
two-face average for each zone was not very different from the heat transfer co- 
efficient for zero angle of attack They concluded that the ovarall heat transfer 
performance of an annular fin is not significantly affected by small departures of 
the angle of attack from zero and since the difference between the heat transfer 
coefficients for the two faces of the fin were found to be significant, measurement 
of the coefficient at only one face may yield erroneous results for the overall heat 
transfer performance 

2.3 Numerical Methods for Solving Navier-Stokes 
Equation 

Different solution techniques have been developed during the past three decades 
The solution techniques for incompressible flows do not have any explicit equation 
for pressure and so due to spatial coupling of pressure and velocity it becomes 
quite difficult to obtain direct numerical solutions For incompressible flow prob- 
lems, pressure does not have the usual thermodynamical meaning Here it is a 
relative variable, which adjusts itself instantaneously for the condition of zero 
mass divergence to be satisfied at all computational cells The primitive variable 
formulation of the incompressible Navier-Stokes equations has following distinc- 
tive feature It is well-known that in an incompressible fluid the speed of the 
sound IS treated to be infinite As a consequence, the pressure field cannot be 
calculated by an explicit time advancement procedure and so it requires at least 
a partially implicit determination which takes into account the coupling between 
the pressure and velocity fields as well as the effects of the velocity boundary 
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conditions This aspect is the most distinctive featuie of the primitive vanable 
foimulation of the incompressible Navier-Stokes equations This pressure field 
can be conveniently computed using the strem-function-vorticity approach for 
two-dimentional flows but it again creates problem for computation of three- 
dimentional flows due to absence of single scalar stream function 

Primitive variable approach can be employed for three-dimentional problems 
which can be followed without encountering non-physical wiggles in pressure dis- 
tribution As a remedy, it has been suggested to employ a different grid for 
each variable Harlow and Welch (1965) have used such a staggered grid for the 
dependent variables in their well known MAC (Marker and Cell) method The 
MAC method of Harlow and Welch is one of the earliest and widely used explicit 
methods for solving the full Navier-Stokes equations In this method, solution of 
velocities are obtained using a two step procedure In the first step called the pre- 
dictor step, the provitional values of velocity components are computed explicitly 
using advection, diffusion and pressure gradient of the earlier time steps This 
explicitly advanced provisional velocity field may not necessarily ensure a diver- 
gence free velocity field Hence, in the second step called the corrector step, the 
pressure and velocity components are corrected so that the velocity field satisfy 
the continuity equation This is done through the solution of a Poisson equation 
of pressure A relative technique, known as pseudo-compressibility method, de- 
veloped by Chorin (1967) involves a simultaneous iteration of the pressure and 
velocity components. Vicelli (1971) has shown that the method due to Chorm 
(1967) and the MAC method are equivalent The original version of the MAC 
method has been modified by Harlow and Amsden, popularly known as Simpli- 
fied MAC (SMAC, 1970) 

Since implicit methods have no such restrictions, they are more attractive Patankar 
and Spalding (1972) have introduced an efficient method known as SIMPLE 
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(Semi-Imphcit Method for Pressure Linked Equation) This method is based 
on a finite- volume discretization of the governing equations on a staggered grid 
In order to improve the convergence involved in the pressure- velocity coupling, 
several variants of SIMPLE algorithm have been developed The SIMPLER al- 
gorithm of Patankar (1981) and the SIMPLEC algorithm of Van Doormaal and 
Raithby (1984) are improvements on SIMPLE Although the changes to incorpo- 
rate SIMPLEC into SIMPLE algorithm are minor, the consequences can be great 
as it eliminates the approximation made in SIMPLE while deriving the pressure- 
velocity corrections 

Application of finite volume methods using non-orthogonal coordinates and col- 
located grids IS reported by Rhie and Chow (1983) and Peric (1985) Ferziger 
and Peric (1999) have shown that collocated arrangement converges faster than 
the staggered variable arrangement and has advantages when extensions such as 
multigrid techniques and non-orthogonal grids are considered Mukhopadhyay 
et al (1993) has developed a numerical method for predicting viscous flows for 
incompressible geometries Integral mass and momentum conservation equations 
are developed and these are discretized into algebric form through numerical 
quadrature The physical domain is divided into a number of non-othogonal 
control volumes which are iso-parametrically mapped on to standard rectangu- 
lar cells Numerical integration for unsteady momentum equation is performed 
over such non-orthogonal cells The results exhibit good accuracy and justify 
the applicability of the algorithm Kobayashi and Pereira (1991) had modified 
the momentum interpolation method suggested by Peric (1985) and named it as 
Pressure-Weighted Interpolation Method-Corrected (PWIMC) In this method, 
the non-orthogonal terms in the momentum equations were solved explicitly, 
whereas in the pressure-corrections they were droped 

It IS clear that subtantial progress has been made for the development of al- 
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gonthms for complex geometries incompressible flow simulations but none of the 
prescriptions is universal Depending on the nature of flow and geometry etc , 
one can always go for the best suited discretization procedure It is very dif- 
ficult to conclude about any algorithm with its universal sense of applicability 
The finite volume algorithm developed by Eswaran and Prakash (1998) for solv- 
ing time dependent three-dimensional full Navier-Stokes equation has following 
advantages It uses collocated grid arrangement, hence all the advantages of col- 
located grid arrangement are gainfully utilized The essence of the algorithm is 
SMAC so it IS suitable for modeling unsteady flows (Kim and Benson, 1992) 
The governing equations are discretized m the physical plane itself without coor- 
dinate transformation This algorithm has been used to solve the flow problem 
and extended to solve energy equation for the present work 

Peric et al (1988) carried out a comparative study of finite volume numerical 
methods with staggered and collocated grids They concluded that the collo- 
cated method possesses no disadvantage relative to the staggered grid version 
The treatment of boundary conditions and the implementation of higher order 
difierencmg schemes are simpler In some cases, the collocated version provides 
faster convergence. It behaves similar to the staggered version under the vari- 
ation of grid and under-relaxation parameters, and provides solutions of equal 
accuracy at the same or lower computing cost 


2.4 How the Current Problem Benefits form Ex- 
isting Literature 

The experiments of Sparrow and Kang (1985) and Apelt et al (1973) reveal the 
effective utilization of the splitter plate for the purpose of enhancement of heat 
transfer while producing less pressure drop Sparrow and Kang (1985) considered 
longitudinal fins at the front, at the rear or both at front and rear The longi- 
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tudinal fin attached at the rear are basically wake splitters The splitter plate 
creates a streamlined extension of the circular tube It brings about enhancemnet 
of heat transfer as well A reduction in the size of wake zone is observed Nar- 
rowing of the wake zone reduces convective heat transfer from the tube surface 
but splitter plate itself generates an extra fin area for conduction Overall, there 
IS an improvement in heat transfer past the circular tube with integral splitter 
plate as compared to the case of flow past a plain circular tube 

The experiments of Sparrow and Chastian (1986) reveal that an annular fin at- 
tached to a tube is not significantly affected by small departure of the angle of 
attack from zero as the two-face averaged heat transfer coefficient does not differ 
from the heat transfer coefficient for zero angle of attack 

The air-cooled condensers of geothermal plants consist of the same fin-tube ar- 
rangements as shown in Figure 1 1 Large fans are used to force the air through 
several rows of these fin-tubes The overall cost of these type of plants is very 
much effected by this type of very large condensers In this present numencal 
study, two three-dimensional numerical models have been formulated to provide 
a better understanding of the heat and flow characteristics The complete Navier- 
Stokes equations together with the energy equation are solved for an element of 
heat exchanger Figure 1 3 is the schemetic representation of the module used for 
the case of a circular tube with an integral splitter plate In this case, the flow 
structure and heat transfer charactenstics for three different lengths (2if, ZH and 
AH) of the splitter plate have been analysed Figure 1 4 shows the module used 
for the case of a circular tube with an attached annular fln The influence of an 
annular fin, attached to the tube at the height of the mid-plane of the channel, 
have been studied 




Chapter 3 


FORMULATION OF THE 
PROBLEM 


3.1 Introduction 

The need to enhance the heat transfer in fin-tube heat exchangers so as to re- 
duce their size and thereby cost, has led us to compute the flow and heat transfer 
characteristics in a horizontal channel having a built-in circular tube with a longi- 
tudinal fin and annular fin. This numerical study can provide information about 
the heat transfer enhancement m fin-tube heat exchangers and effect of various 
flow parameters 

3.2 Statement of the Problems 

3.2.1 Splitter plate 

The computational domain for the present problem is shown in Figure 1 3 Com- 
putations are performed for a circular tube of diameter JD = 4 OH, H being the 
height of the domain parallel to axis of the tube The domam length along flow 
direction is Li = 20if and the width in transverse direction is L 2 = 11 25R 
The center of the tube is at a distance Ic = 0 4R from the inlet For a par- 
ticular flow Reynolds number, different fin lengths are considered (I = 2H, 3H 
and 4H) for flow and heat transfer study. Again, for a particular fin length. 
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computations are performed to investigate the effect of variation of Reynolds 
number(J?eH=500, 1000 and 1500) on the behavior of heat transfer Air has been 
considered as the working fluid, hence the Prandtl number is taken as 0 7 The 
tube and the longitudinal fin (splitter plate) are considered to be isothermal sur- 
faces 

3.2.2 Annular Fin 

The computational domain is shown in figure 1 4 Computations are performed 
for a circular tube of diameter D = 4 QH, H being the height of the domain 
parallel to axis of the tube The domain length along flow direction is Li = 25 H 
and the width in transverse direction is L 2 = 12 dH The center of the tube is 
at a distance Ic = 8 from the inlet The outer diameter of the annular fin is 
£>1=1 5£> and Reynolds number is i?ei3-=1000 for flow and heat transfer study 
Air has been considered as the working fluid, hence the Prandtl number is taken 
as 0 7 


3.3 Governing Equations 


The three-dimensional continuity, Navier-Stokes and Energy equations for lami- 


nar flow of an incompressible fluid having uniform density are 


du^ 

dxi 


= 0 


( 31 ) 


dut du^ _ 1 dp d'^Ut 

dt dxj p dxi ^ dxjdxj 

_ d^T 

dt '^^^dxj ^dxjdxy 


(3 2) 
(3 3) 


3.4 Boundary Conditions 


• Top and bottom walls 

These are considered as no-slip walls for x and y-components of velocities 
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u = V = w = 0 and = 0 
r = {T^ represents wall temperature} 

• Side walls 

dy ~ aj/ “ 

1 = 0 and f = 0 

• Channel inlet 

u = Uoo, V = w = 0 and ^ = 0) 

T = Too 

\ 

• Channel exit 

The mass flux through the outlet boundaries is found by by means of a con- 
tmuative outflow condition (Orlanski (1976)), which allows changes inside 
the flow field to be transmitted outward, but not vice versa. This boundary 
condition uses an upwmded form of the following equation for each velocity 
component and temperature 
^ + Uav^ = 0 (where (j) represents u,v,w ovT) 

The pressure is also specified at the outflow 
P=Poo 

• Obstacles (Circular tube. Splitter Plate and Annular Fin) 

u = V = w = 0, ^ = 0 (where n signifies the normal direction to the 

surface) 

and T = Tyj 




Chapter 4 


GRID GENERATION AND 
SOLUTION TECHNIQUE 

4.1 Introduction 

The numerical solution of the governing equations in the prescribed computa- 
tional domain with given initial and boundary conditions involves two major 
discretization steps First of which is the discretization of the domain with the 
help of a numerical grid The governing equations are discretized on the basis of 
this numerical grid Various ways of discretization can be applied to the above 
set of partial differential equations Among these the most common methods are 
those using finite differences, finite element and finite volume In the present 
work a finite volume scheme is chosen for spatial discretization of the governing 
equations As a result, a conservative formulation of the discretization is possible 
This automatically assures the satisfaction of the global balance of the conserved 
quantities, independent of the coarseness of the numerical grid. 



Figure 4 1 The Grid System for Splitter Plate 



Figure 4 2 The Grid System for Annular Fin 
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4.2 The Grid 

Figure 4 1 and 4 2 are the schematic presentations of the two-dimensional grids 
used respectively in the cases of a circular tube with an integral splitter plate 
and a circular tube with an annular fin A block partitioning method is used to 
generate the grid The two-dimensional grid-mesh is stacked m the third direc- 
tion to obtain the three-dimensional grid-mesh The grid generation technique 
adopted for body-fitting type is described in brief in the following sections Three 
categories of physical domains may exist viz (1) a simply-connected domain (2) 
a double-connected domain and (3) a multiply-connected domain Here simply- 
connected domain is used in the body-fitting region 

4.2.1 Grid Generation Technique 

The two-dimensional grid has been generated using tranfinite interpolation and 
dififerential equation method Transfinite interpolation uses linear interpolation 
scheme to compute the interior points by using the coordinates of the boundaries 
The grid produced by this method is not enough smooth because of propagation 
of discontinuities from the boundary to interior points of the domain That is 
why the mesh obtained by algebraic mapping is further improved by the use 
of Partial Differential Equations (PDE) technique in which a system of PDE is 
solved for obtaining the location of the grid points m the physical space But in 
the computational space the transformed grid is of rectangular shape with uni- 
form spacing We can use any of the three types of partial differential equations, 
namely elliptic, parabolic and hyperbolic Elliptic equations are prefered for 
closed geometries, hyperbolic equations are used for domain where outer bound- 
ary is not prescribed and parabolic equations are used where boundary is closed 
on one side (prescribed) and open on the other side The present grid is generated 
using of elliptic partial differential equations In any gnd generation technique 
geometric information is fetched from the boundanes Hence the steady state 
boundary-valued nature of elliptic ^uations make them most favourite and so 
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they are widely used m grid generation A Laplace equation or a Poisson equation 
with Dirchlet boundary condition can be used for this purpose Standard Poisson 
equation is of the form 

V2^=P 

Where P and Q are known as control functions If P=Q=0 then, the equations 
get modified to Laplace form as 

V2^ = 0 
V^ri = 0 

The above equations can be solved by finite difference technique to get the loca- 
tion of the interior points from the boundary The initial guess needed to solve 
the Laplace equation can be generated using algebraic mapping This Laplacian 
operator can provide quite smooth grid The grid lines remain equally spaced 
if the boundary curvature is absent but near the curved boundaries they show 
a tendency to concentrate This drawback can be controlled with appropriate 
selection of control functions P and Q Also control functions can be objectively 
employed for enforcing orthogonality over the selected edges. That may simplify 
implementation of boundary conditions, but use of control functions remains as- 
sociated with additional complexity in the transformed domam 

4.2.2 Computational Procedure for Grid Generation 

The algorithm used for the generation of the present grid can be outlined as 
follows 

1 Geometrical data is used as input. 

2 The computational grid is defined on ^ plane based on the number of gnd 
points 

3 Boundary points along the edges of the physical domain are defined 

4 Algebraic gnd is used to generate the initial guess by hnear interpolation 
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technique, which is used by the Laplacian operator 

5 Using four-point formula, guess is further improved a little 

6 Elliptic grid generation system constitutes of the following equations derived 
from the Laplace operator 


ax^^ - 

- 2bx^ri + cXnn = 0 

(41) 


- 2by^fi cyfjTj = 0 

(4 2) 


where 


a = + y2 

b = x^yr, ■+- Xr,y^ 
c-x(^ + y^^ 

Here subscripts indicate partial derivatives of cartesian variables with respect to 
curvilinear variables 

These coefficients are determined using finite difference approximations Their 
values in these expressions are provided by the initial distribution for the first 
iteration, and subsequently from the previous iteration, i e the computation of 
coefficients lag by one iteration level The iterative solution continues until a 
specified convergence criterion is meet 

Because of discretized approximation of differential equations, the truncation 
error in numerical analysis is much dependent on the nature of gnd We target 
for a grid, which gives minimum truncation error. 

The parameters that affect the quality of the grid are 

1 Transformation Jacobian 

2 Skewness 

3 Aspect Ratio 

4 Adjacent cell ratio 
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A desirable grid should have nonzero transformation Jacobian almost at each 
location of the grid cells It should be of similar order of magnitude with no 
quite small values appearing anywhere It should possess least possible skewness, 
aspect ratio close to 6 and adjacent cell ratio close to one 
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4.3 Solution Technique 

The present work uses a finite volume scheme based on the conservative formu- 
lation, assuring the gobal balance of the conserved quantities, independent of the 
coarseness of the numerical grid The conservation equations are dicretized by 
employing the finite volume approach of Eswaran and Prakash (1998) and Prab- 
hakar et al (2002) The solution domain is divided into a number of contiguous 
(finite) control volumes (CV) The control volumes are defined by the coordi- 
nates of their vertices, which are assumed to be connected by straight lines The 
coordinates of the control volume vertices are calculated by the grid generation 
procedure A Collocated Grid arrangement is employed and all the dependent 
variables u,v,w,p and T are defined at the same location - the centroid of the 
control volume (Figure 4 3) The symbols E, W, N, S, T and B indicate the 
SIX neighbouring control volume centers for the east, west, north, south, top and 
bottom neighbours respectively 

Discretization of General Equation 

(a) Rate of Change 

The value of the dependent variable <f> at the centroid of the control volume (the 
geometric center) represents an average over the CV as a whole Thus 


dt 



ipw)r^ - jp^v); 

At 


pV 


At 


where V is the volume of the cell 


(4 3) 


(b) Convection Fluxes 

The surface integral over convection flux of variable (j) can be approximated in 
the following form 
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pu(j> dS p(f>j (u S)j = Y (4 4) 

J 3 

where <^j is the value of (j) at the center of face j Thus 

pil(j> dS Fe(j)e + Fyj(j>ju + Fn<f>n + Fs(j)s + Ft(j>t + (4 5) 

where <^e is the interpolated value of the variable (j) at the east face center and 
so on This can be evaluated by using an interpolation scheme which is based on 
blending of fluxes due to upwind and central difference schames 

(c) Diffusion Fluxes 

Ferziger and Peric (1999) developed an expression for the diffusion fluxes for the 
finite volume approach Eswaran and Prakash (1998) also explained the means 
of expressing diffusion fluxes on different faces of arbitrary-shaped differential 
control volumes Following Eswaran and Prakash (1998), the diffusion flux of a 
vanable (j> through the cell faces is evaluated in the following manner as follows 



/’r*v,*ds« = 

'' ^ j=e,w,n,s,t,b J 


For any face we can write 


(4 6) 


Sj = + o.2X^ CKsn^ (4 7) 

where n^, and are any three linearly independent (not necessarily orthog- 
onal) unit vectors Therefore, 


(4 8) 


V^.Sj = V(f> (ain^ -f 0 : 211 ^ + 0 : 311 ^) 

= Q-iV^ + a2V(f> + azV(j> n® 
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The procedure for calculating ai, and was described by Eswaran and 
Prakash (1998) and Prabhakar et al (2002) The diffusion flux is made up 
of two distinct parts - noimal derivative diffusion flux and cross derivative dif- 
fusion flux The second part arises from the non-orthogonahty of the grid The 
example of the east face is taken to illustrate the diffusion model (Figure 4 4) 
Given the edge center values 4>te, (pbe-, 4>se, <t>ne we can get the normal diffusion 
term and the cross diffusion term and Finally, the diffusion 

flux is computed by 




<j>E 




4*p I ^te 
+ CX2 


Ax^ 


+ 0:3 




Arc® 


(4 9 ) 


Pressure and Velocity Coupling 

To obtain the velocity and pressure field satisfying the mass and momentum 
conservation laws, a procedure akin to the Marker and Cell method is used 
This method offers an efficient and easy way of pressure- velocity coupling It is 
basically a semi-explicit method The momentum equations are discretized in an 
explicit manner with the exception of the pressure gradient terms that are treated 
implicitly, and the continuity equations are also enforced implicitly This can be 
expressed mathematically by the following two discretized equations 


pV 


u' 


n+l 


U'i 




At 


( 410 ) 
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The momentum equations are solved using the guessed values of the velocity and 
pressure field The provisional velocity components u* are calculated from the 
following equation 

u;. = u”p-^(F|+F^) + ^S. (412) 

where u" is the value of the velocity at earlier iteration and is the pressure 
term 

This provisional velocity in general will not satisfy the continuity equation 
The continuity Equation 4 11 m another form reads as follows 

^(F;+i/) = 0 (413) 

3 

3 3 

where the right side of the Equation 4 14 is the uncorrected mass flux obtained 
from the provisional velocities and the left hand side represents mass flux correc- 
tion To evaluate the terms on the right hand side of Equation 4 14, it is necessary 
that the variables (velocity and pressure) are known at each cell faces Owing 
to the non-staggered arrangement, the linear interpolation between the adjacent 
cell centers leads to a checker board pressure field as shown by Rhie and Chow 
(1983) This problem can be circumvented by using the concept of momentum 
interpolation The pressure correction equation finally becomes 

AtVpj S = pVySj err F'^ = -AtVp' (4 15) 

After the soluion of the pressure correction equation, the nodal velocity, mass 
fluxes and pressure are updated Using the converged velocity, the energy equa- 
tion IS discretized and solved by a SOR tedbnique to find the temperature distri- 
bution. 



Chapter 5 

Results and Discussion 


5.1 Introduction 

This chapter is divided into two parts Flow and heat transfer characteristics for 
the case of a circular tube with the splitter plate are discussed in the first part 
The second part describes the flow and heat transfer for the case of a circular 
tube with integral annular fin 

A 61 X 49 X 21 grid-mesh is used for the case of a tube with a splitter plate 
Computations have been carried out for three different lengths of the splitter 
plate (L = 2H, 3H and iH) and three different Reynolds numbers (500, 1000 
and 1500) The length of the splitter plate for all other parametric variation is 
considered to be T = 2H except for the case when the comparison of heat transfer 
performance is made for different lengths of the splitter plate 

A 60 X 41 X 17 grid-mesh is used for the annular fin case. Computation has 
been earned out for annular fin of outer diameter Di = I bD and a Reynolds 
number of 1000 

The divergence-free criterion is satisfied using an upper bound of 10~^. Span- 
averaged Nusselt number is calculated on the basis of bulk-mean temperature 
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Air IS considered as the working fluid and hence Pr = 0 7 


5.2 Splitter Plate 

The use of splitter plates to add surface area can bring about substantial increases 
in the tube-bank heat transfer The presence of the splitter plate also influences 
the pressure drop characteristics of the tube bank 

5.2.1 Flow Characteristics 

Figure 5 1 shows the particle paths of an instantaneous flow fleld at the horizontal 
midplane for flow past a circular tube and the circular tube with an integral split- 
ter plate In Figure 5 1(a), the flow past a circular tube displays vortex shedding 
It IS evident from the asymmetry observed in the wake structure It is already 
established in literature that the vortex shedding for flow past a circular tube, 
placed in an infinite medium, begms at a low Reynolds number {Re ~ 90) The 
Reynolds number in the present study {Ren^lOOO) is quite large and alternate 
shedding of vortices is an ubiquity The wake splitter has been frequently used as 
a passive means of controlhng vortex formation and sheddmg in the near wake of 
a circular tube The flow field at the horizontal mid-plane for the case of flow past 
a circular tube with integral splitter plate (Fig 5 1(b)) illustrates a symmetnc 
behavior even though it is an instantaneous flow field For the plain circular tube 
(without a splitter plate) the shear layers exhibit transverse oscillations in-phase 
with the Karman shedding cycle. The attachment of even a very short splitter 
plate IS sufficient to di mmis h the amplitude of the oscillations The amount of 
vorticity shed from the circular tube m each cycle is not affected by the splitter 
plate since the plate does not alter the vorticity produced in the boundary layer 
or at the point of separation 


Figure 5 2 shows the streamlines on the horizontal midplane of the channel for 
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two different lengths of the splitter plate (L = SH and L = 4H) In both the 
cases the flow fields illustrate a good degree of symmetry 

A plain circular tube has its vortex formation region well forward as m Figure 
5 1 (a) Splitter plate of length 2H causes the vortex to form clear of the plate 
as m Figure 5 1 (b) As the length of the plate increases, vortices move towards 
the plate, roll over the edge as in Figure 5 2 The position of vortex formation 
moves downstream as L increases from 2H to 4H, being near the plate trailing 
edge in each case Thereafter the formation position remains at essentially the 
same location with respect to the circular tube The formation region at L = AH 
IS slightly downstream of all the others 

Figure 5 3 shows the limiting streamline plots m the region close to the plate 
for a plain circular tube and the circular tube with splitter plate of length =2H 
for a Reynolds number of 1000. The two plots are quite similar except for the 
zone where the splitter plate is attached A saddle point of separation and a 
horseshoe vortex system are observed from the figure The incoming flow does 
not separate in the traditional sense but reaches a stagnation or saddle point of 
separation and goes around the body The noddle point of attachment (marked 
C) and the separation lines which from circular arcs across the tube are seen in 
the figure The flow above the lower wall hits the front of the tube A signif- 
icant part of it moves downward and creates a region of reversed flow m front 
of the stagnation line On each side of the tube, there is a region of converging 
streamlines (marked G) These are the traces of horseshoe vortices The mean 
shear within the approaching boundary layer is skewed or deflected by the trans- 
verse pressure gradient Accordingly, the boundary layer separates and rolls up to 
form a spanwise vortex The horseshoe vortex has high influences on the overall 
convective heat transfer There is a wake stagnation point (marked F) further 
downstream of the body 
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The structure of the surface streamlines or limiting streamlines on the surfaces 
of both circular tube and splitter plate has been investigated These streamlines 
consist of only tangential components of velocity near the surfaces of the tube and 
the splitter plate Figure 5 4 and 5 5 show such streamlines The time-averaged 
flow fleld is symmetric and the limiting streamlines are symmetrically distnbuted 
about the midplane of the channel The surface streamlines reveal bifurcation 
The bifurcation could be either a positive bifurcation or a negative bifurcation 
When one streamline divides into two or more streamlines, it is called a positive 
bifurcation If two or more streamlines combine to form a single streamlme, it is 
called a negative bifurcation Figures 5 4 and 5 5 reveal the positive bifurcation 
at the forward stagnation line of the tube surface (0 = 0 or 9 = 2%) It 
IS also observed from the Figures 5 4 and 5 5, that for a plain circular tube the 
line of reattachment is at the 0 = tt and for the splitter plate it is at the tip of 
the plate As the length of the splitter plate increases, the streamlines on either 
side of the line of reattachment are straightened Figure 5 4(a) shows two other 
lines at 0.6l7r < 0 < 0 SItt and 1 IQtt < 0 < 1 SOtt respectively Near these hues, 
the tangential components of the shear stress vector vanish Hornung and Perry 
0 called this a negative stream surface bifurcation These lines can be termed 
separation lines, along which the boundary layer separates from the tube surface 
The separation angles (on the tube surface) for the cases with splitter plates (of 
varying lengths) remain almost same when compared to the separation angles 
for the case without a splitter plate. However, this is not directly evident from 
Figures 5 4(a) and 5 5 since the lengths of the splitter plates have been added 
to the scales for the abscissa. A closer observation reveals that the presence of 
splitter plates is slightly advances the flow separation 


Figure 5 6 shows the distnbution of pressure coefficient Cp — j where 

p IS the local pressure, is the stagnation pressure and Poo is the atmospheric 
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pressure for the case of a circular tube in uniform flow with splitter plate of 
diflterent lengths at Ren = 1000 The profiles indicate that the pressure in the 
separated region, adjacent to the circular tube is highly sensative to the variations 
in length of the splitter plate However the flow separation point is not altered 
by the presence of the splitter plate These findings agree well with the trends of 
variation in the earlier work of Apelt et al (1973) The reason for an increase in 
the base pressure due to presence of the splitter plate appears to be as follows It 
was mentioned earlier that the splitter plate advances the flow separation slightly 
However this early separation does not cause any adverse effect in the wake zone 
The small rise in base pressure is expected to decrease the total drag 
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Figure 5.2 


Streamlines on the horizontal midplane of the channel for 

(a) circular tube with splitter plate of Lengm=3H 

(b) circular tub© with splitter plate of Lengui— 4H 
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Figure 5.3 Streamline very near to bottom plate of the channel for 

spBtler Pbte of length=2H 
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Figure 5.4 Limiting streamline on the surface of (a) circular tube 
and (b) circular tube with splitter plate of L=2H 
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Figure 5.5 Limiting streamline on the surface of 

(a) circular tube with splitter plate of length=3H 

(b) circular tube with splitter plate of length=4H 
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Figure 5.6 Comparison of distribution of Cp around the tube surface 
for four different cases ( L=0, L=2H, L=3H and L=4H ) 
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5.2.2 Heat Transfer Characteristics 


Temperature Distribution 

Figure 5 7(a) shows the temperature distribution of the cold fluid at the loca- 
tions close to the tube surface for the plain circular tube Figure 5 7(b) shows 
the same for the case of a circular tube with a splitter plate of length =2H It 
is observed from the Figure 5 7(a) the mmimum temperature is at the forward 
stagnation point The temperature increases as the flow moves towards the point 
of reattachment but there is a slight decrease m temperature just before this reat- 
tachment point It could be explained in the following way. The incoming cooler 
fluid strikes the stagnation line and removes heat from the stagnation region Af- 
ter staking the tube, the fluid gams some thermal energy and this energy in the 
fluid increases by taking heat from the tube surface by convection as it progresses 
embracing the tube surface This flgure also shows that the temperature is higher 
near the bottom and top plates as compared to the middle portion of the channel. 
Now at the splitter plate part it is observed that as the fluid progresses towards 
the tip of the plate, its temperature decreases and at the tip it becomes minimum 
From the Figures 5 7(b) and 5 8 it is observed that as the length of the splitter 
plate IS increased, the low temperature zone dominates in the z direction. The 
effect of splitter plate on the temperature field becomes more dominant as the 
length of the plate increases The splitter plate of length 2H and 3H do not have 
much influence near the bottom and top plates 

In these two figures, splitter plate of different lengths show symmetry with re- 
spect to z=H/2 through its instantaneous field The reason could be splitter plate 
imposes symmetry by streamhmng the flow 
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Nusselt Number Distribution 

Figure 5 9 shows the contour map of Iso-Nusselt number on the bottom plate 
Before starting the discussion on this plot, there are some issues which are to be 
addressed In experiments, the tube-wall and the fin do not have constant tem- 
perature because both convective and conductive heat transfer take place from 
the tube-wall and the fin But in the present computation, it is assumed that the 
tube-wall and the fin, are at the constant temperatue It is also to be noticed that 
in the present numerical simulation, the side walls are considered as free-slip walls 
whereas m the experiments, the side walls are of no-slip type Figure 5.9 depicts 
that at the leading edge of the bottom plate, Nusselt number has its maximum 
value and it decreases gradually till some specified position before the tube At 
the leading edge the cooler fluid comes m contact with the hot solid-wall for the 
first time hence the temperature gradient becomes maximum entailing very high 
value of the Nusselt number The inlet condition is one of uniform velocity as well 
as uniform inlet temperature Therefore, both the velocity and thermal boundary 
layers are developing at the channel leading edge The region that follows the 
leading edge of the channel is the combined entrance region There is a sudden in- 
crease m Nusselt number in front of the tube. The reason for the abrupt increase 
in Nusselt number can be attributed to the formation of horseshoe vortices The 
spiraling motion of the horseshoe vortices brings about a better mixing and the 
heat transfer in this region is enhanced significantly The Nusselt number is low 
in the wake region The separated dead water zone where fluid recirculates with 
low velocity causes poor heat transfer Figure 5 9(b) and 5 10 show the contour 
plots of iso-Nusselt number distibution on the bottom plate for the splitter plate 
of length (L = 2H, SHandAH) The Nusselt number distribution pattern from 
the leading edge to the location of the tube in all the cases is sinular to that of 
the plain circular tube The heat transfer enhancement due to the attachment 
of splitter plate is evident m the zone near the tip of the plate As the length of 
the splitter plate increases, the size of the zone of high heat transfer in the wake 
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increases 

The presence of splitter plate causes a general decrease in the local heat transfer 
rate over the surface of the circular tube The decrease in local heat trasnfer is 
caused since the interaction of the vortices on either side of the splitter plate is 
ceased However, the splitter plate provides with additional heat transfer sur- 
face and overall heat transfer is enhanced A parameter could be formulated to 
estimate the gam in heat transfer due to the splitter plate in the following way 

^ L. (5 1) 

•^0 

Here, Aq is the surface area of the tube and A is the combined surface area of 
the tube and the splitter plate 

Table I Values of Q/Qo for splitter of three different lengths 


Sphtter plate length (L) 

A/Ao 

Q/Qo 

2H 

13183 

1 196 

3H 

14775 

1277 

4H 

1 6366 

1369 


The computaions are performed for Re = 500, 1000 and 1500 and the length of 
the splitter plate L = 2H,3Hand4:H. The Nusselt number {Nuavg) for different 
lengths of the splitter plate is calculated for a Reynolds number of 1000 The 
Q/Qo ratios deduced from these Nusselt numbers axe shown in Table I The per- 
tinent parameters of the Table axe L, A/ Ao and Q/Qo As to symbols, Q is the 
heat transfer from Hie combined surface area of the splitter plate and the circular 



tube and Qq is the heat transfer from the circular tube alone It is clear from 
the table that Q/Qo > 1 when A/Aq > 1 For the case of = 1 60, (Q/Qo 
) could be compared with the results of Sparrow and Kang(1985) In our com- 
putation we observe Q/Qo = 1 369 which is less than the result of Sparrow and 
Kang (1985) by 7% The reason could be as the following The tube length in the 
experiments is quite large so that a nominally two-dimentional situation evolves 
and transport mechanism is free from the influence of top and the bottom walls 
In the present computation, the top and the bottom walls play a significant role 
The gap between the top and the bottom wall is quite small and the length of 
the tube is restricted by the gap 


Table II Variation of Nuavg with Reynolds number for splitter plate of 

length=2H 


Reynolds number (Ren) 

^ ^avg 

500 

4 535 

1000 

5.23 

1500 

6 038 


The Table II shows the variation of Nuavg with Reynolds number. As the 
Reynolds number increases the Nuavg also increases Increase m Reynolds num- 
ber means increase in inlet velocity It is understood that the splitter plate is more 
effective at higher Reynolds numbers 


The splitter plate creates a streamlined extension of the circular tube It brings 
about enhancemnet of heat transfer as well A reduction m the size of wake zone 
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IS observed Narrowing of the wake zone reduces convective heat transfer from 
the tube surface but splitter plate itself generates an extra fin area for conduc- 
tion Overall, there is an improvement in heat transfer past the circular tube 
with integral splitter plate as compared to the case of flow past a plain circular 
tube. The heat transfer enhancement obtained by finning can be compared with 
that which could be obtained by increasmg the diameter of unfinned tubes The 
computation has been carried out for equal transfer surface area of unfinned tube 
having increased diameter and finned tube at a particular Reynolds number of 
Rep =4000 The area ratio considered to investigate the performance of finned 
tube banks is A/A.o=l 3183 At that area ratio of 1 3183 and for Rep=4000, 
the value of Nuavg is 5 23 for the case of a circular tube with splitter plate of 
length 2H and the value of Nu^vg is 4 88 for the case of a plain circular tube of 
increased diameter So finning gives better result for heat transfer as compared 
to the increased diameter circular tube 

The Figure 5 11 shows the comparison of span-averaged pressure variation be- 
tween finnned tube and unfinned tube of same surface area This comparison 
is presented for a fixed value of Reynolds number (i?ep=4000) The figure indi- 
cates that the pressure penalty is very high m unfinned circular tube compaired to 
finned tube Here increased diameter circular tube provides large sized obstracle 
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Circumferential length 

Figure 5.7 Temperature distribution of the cold fluid near the surface of 

(a) circular tube 

(b) circular tube with splitter plate of length=2H 



i 


Nu 




29 2717 
27 3519 
25 4321 
23 5123 
21 5925 
19 6727 
17 7529 
15 8331 
13 9133 
11 9934 
10 0736 
8 15383 
6 23402 
4 31422 
1 99325 



Nu 




29 4451 
27 4851 
25 5251 
23 5651 
21 6052 
19 6452 
17 6852 
15 7252 
13 7652 
11.8052 
9 84526 
7 88527 
5 92529 
3.96531 
2 00533 


Figure 5.9 Iso-Nusselt number distributatlon at the bottom plate for 
fs) circulsr tub© 

(b) circular tube with a splitter plate of length=2H 
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Figure 5.1 0 Iso-Nusselt no. distribution at the bottom plate for 

(a) splitter plate of length=3H 

(b) splitter plate of length=4H 
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Figure 5.1 1 Comparison of span avaraged pressure distribution 
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5.2.3 Specific Studies for the Pertinent Engineering Pa- 
rameters 

In order to reveal the impotance of this problem on engmermg applications, the 
following parametric study has been undertaken The coefficient of form drag 
(Cdf) and coefficient of Viscous Drag (Cdv) have been calculated for the varying 
length of the splitter plate The averaged Nusselt number (Nu^avg) oJily on the 
tube surface has been computed for different lengths of the splitter plate Also 
the global span averaged Nusselt number (Nu^avg) at the bottom plate has also 
been calculated for the varying geometrical parameters Table III summarizes the 
study of drag coefficients as the function of length of the splitter plate Likewise 
Table IV depicts the variation of tube surface based averaged Nusselt number 
and global Nusselt number on the bottom plate vis-a-vis the length of the splitter 
plate 


Table III Variation of Cdf Cdv 


Length of the Splitter Plate 

Cdf 

Cdv 

Total Drag 

0 

1234 

0 0453 

1 2793 

2H 

0 9911 

0 0433 

1 0344 

3H 

0 9856 

0 0423 

1 0279 

4H 

0 983 

0 0411 

10241 


It can be seen from Table III that, with the increase in length of the splitter plate, 
the value of Cdf and Cdv axe decreasing The value of the Cdv is very less for 

each cases 
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Table IV Variation of Nu^avg and Nu^avg with splitter plate length 


Splitter Plate Length 

•ATt/ Q/ug 

avg 

0 

5 765 

6 388 

2H 

5 603 

6 417 

3H 

5 637 

6 441 

4H 

5 687 

6 462 


Table IV demonstrates, as the length of the splitter plate increases, the averaged 
Nusselt number based on tube surface first decrease and then increases It also 
reveals that the global Nusselt number at the bottom plate increases very slightly 
with the increase in length of the splitter plate 
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5.3 Annular Fin 

Annular fins are of great practical importance for their use in heat exchangers 
For a given fin weight, the fin can dissipate various quantities of heat, depending 
on its shape and geometry 

5.3.1 Flow Characteristics 

Figure 5 12 shows the streamlines on a horizontal plane which is at the midway 
between the annular fin and the channel wall It shows two symmetrical vortices 
in the near wake of the circular tube This plot is similar to the streamline plot on 
a horizontal plane mid-plane of a plain circular tube placed in a channel These 
two vortices are marked as the Karman vortices The flow separation in the wake 
zone IS clearly visible 

Figure 5 13 shows the limiting streamlines the bottom plate for a circular tube 
with annular fln A saddle point of separation and a horseshoe vortex system are 
observed from the figure The incoming flow does not separate m the traditional 
sense but reaches a stagnation or saddle point of separation and goes around the 
body The noddle point of attachment (marked C) and the separation lines which 
from circular arcs across the tube are seen in the figure The flow above the lower 
wall hits the front of the tube A significant part of it moves downward and 
creates a region of reversed flow m front of the stagnation line On each side of 
the tube, one finds a region of converging streamlines (marked G) These are the 
traces of horseshoe vortices The mean shear within the approaching boundary 
layer is skewed or deflected by the transverse pressure gradient Accordingly, the 
boundary layer separates and rolls up to form a spanwise vortex The horseshoe 
vortex has high influences on the overall convective heat transfer There is a wake 
stagnation point (marked F) further downstream of the body 

Figure 5 14 shows the streamlines plot on the vertical mid-plane The figure 
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leveals symmetry about the mid-channel-height &i z = H/2 The figure also 
shows the the footprints of the horseshoe vortices (marked as Si and S^) As the 
flow progresses towards the annular fin and strikes it, some amount of fluid moves 
up and some amount of fluid moves down At the tip of the annular fin, the flow 
reattachment is clearly visible A pair of symmetrical vortices is observed in this 
figure The formation of these two vortices is influenced by the pressure difference 
between the tip of the annular fin and plates 

Figure 5 15 represents the stremlines on the cross strem plane at a location x — lc 
{Ic is distance of the center of the tube from the inlet) 

Figure 5 16 represents the distribution of pressure coefficientC'p= > where 

p is the local pressure, Ps, the stagnation pressure, Poo, the atmospheric pressure 
around a circular tube wall One of the two curves is for Cp around the circular 
tube where annular fin is not attached and another one is for the Cp around the 
a nnular fin The distibution of the Cp around the circular tube surface is same 
as that for the plain circular (without any attachment) tube case It is observed 
from the figure that there is an incomplete recovery of pressure This is due to 
the recirculation behind the circular tube. The point of separation is clear from 
the figure 



Re=1000 Pr=0 70 D/B=0.32 




Re=1000 Pr=0 7 D/B=0.32 



Figure 5.13 Streamline very near 



mid-plane 







Figure 5.16 Distribution of pressure coefficient around the tube with annular fin at 
different positions 
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5.3.2 Heat Transfer Characteristics 

In this section the heat transfer characteristics for the annular fin is discussed 

Temperature Distribution 

Figure 5 17 shows the temperature distribution on the vertical midplane of the 
channel This figure reveals the development of thermal boundary layer along the 
top and bottom plate The contour plot is symmetric with respect to z = JEr/2 
The effect of the annular fin on the temperature contours is observed in this fig- 
ure 

Figure 5 18 represents the temperature contours at the cross strem plane at a 
location X = Ic This figure shows how the plates and the annular fin influences 
the temperature of the fluid 


Nusselt Number Distribution 

Figure 5 19 Represents the contour map of local Nusselt number very close to the 
bottom plate and the top plate These two plots are identical. There are some 
issues which are to be addressed before hand In practice, the tube, channel walls 
and the fin do not have constant temperature because both convective and con- 
ductive heat transfer take place from the tube, channel walls and the fin In the 
present computation, it is assumed that all (tube, channel walls and the annular 
fin) are at the constant temperatue At the leading edge of the bottom plate, 
Nusselt number has its maximum value and it decreases gradually upto some 
location infront of the tube At the inlet, the cooler flmd comes in contact with 
the hot solid wall for the first time and the temperature difference is maximum. 
The Nusselt number reveals highest value herein The thermal boundary layer 
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developes downstream The formation of thermal boundary layer causes decrease 
m heat transfer rate and the Nusselt number decreases There is a sudden in- 
crease in Nusselt number in front of the tube The reason for the sudden increase 
in Nusselt number can be attributed to the formation of horseshoe vortices The 
spiraling motion of the horseshoe vortices brings about a better miymg and the 
heat transfer m this region is enhanced significantly 

The local Nusselt number distribution on the top and bottom side of the an- 
nular fin IS shown in Figure 5 20 The Nusselt number distribution pattern at 
the top surface of the annular fin is same as that for the bottom surface Here, 
at the extreme left, the relatively colder fluid comes m contact with the heated 
fin surface The heat transfer is maximum here As the fluid moves above and 
below around the fln surface, its temperature increases and Nusselt number on 
the surfaces decreases 
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Figure 5.17 Temperature distribution at the plane(Z-X) cut through 
the centre of the tube 
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Figure 5.1 8 Temperature distribution at the vertical plane (Y-Z) cut through 
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Nu 



24 4906 
22 9235 
21 3565 
19 7895 
18 2225 
16 6555 
15 0885 
13 5214 
11 9544 
10 3874 
8 82039 
7 25337 
5 68635 
411934 
2 55232 



Nu 



24 5031 
22 9353 
21 3675 
19 7997 
18 2319 
16 6641 
15 0963 
13 5284 
11 9606 
10 3928 
8 82504 
7 25724 
5 68943 
4.12163 
2 55383 


Figure 5 


.19 Iso-Nusselt number distribution at the 
(a) bottom plate &(b) top plate 
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Figure 5.20 Nusselt number variation of the Annular F 

(a) atthetf)pewaoe , 

(b) at the bottorh surface - : 
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5,4 Model Validation 

The model validation was performed through comparison with the experimental 
results of Eckert and Soehngen (1952) and numerical results of Chun and Boehm 
(1989) for a Reynolds number of 218 Figure 5 21 shows the comparison of local 
Nusselt numbers along the circumference of the tube Present computation com- 
pares favourably with the experiments of Eckert and Soehngen (1952) 




Figure 5.21 Comparison of local Nusseit number distribution on the tube surface 



Chapter 6 


Conclusions and Scope for Future 
Work 

6.1 Conclusion 

A three-dimensional computational study of forced convection heat transfer has 
been carried out to determine the flow and heat transfer characteristics The full 
Navier-Stokes equations together with the energy equation are expressed in mte- 
gral form A semi-explicit control volume formulation of Eswaran and Prakash 
(1998) is used to solve the governing equations This study focuses on the detailed 
analysis of the flow and heat transfer characteristics in a rectangular channel with 
a built-in circular tube havmg a longitudmal fin, i e , Splitter Plate and a built-in 
circular tube having a circumferential fin,i e , Annular Fm 
Particular attention has been paid to the variation of coefficient of pressure, co- 
efficient of drag, vortex structure, limiting streamlmes and heat transfer with 
splitter plate length This study has shown that the characteristics of the wake 
downstream of a circular tube can be considerably altered by placing a splitter 
plate on the wake centerline downstream of the circular tube Flow visualization 
indicated that a splitter plate reduced the level of 3-dimensionality in the forma- 
tion region by stabilizing the transverse ’flapping’ of the shear layers The use 
of splitter plate streamlmes the flow Thus the vortices are pushed downstream, 
which decreases the heat transfer But again the presmce of sphtter plate in- 




( 'f. 


If 


6 2 SCOPE FOR FUTURE WORK 


70 


creases the heat transfer area substatially by complementing for conduction heat 
transfer This results m an overall enhancement m heat transfer The splitter 
plate being a slender body m the shown configuration, tends to decrease even 
the pressure loss penalty A circular tube with an annular fin causes a significant 
change in the temperature field of the channel The presence of annular fin in- 
creases the heat transfer area The Nusselt number distribution on the top and 
the bottom surfaces of the annular fin shows the pattern which is quite similar 
to that of bottom and top wall of the channel The pressure distnbution around 
the annular fin is quite similar to that of a circular tube 


6.2 Scope for Future Work 

The results are presented for the case of splitter plate of lengths {L = 2i?, 3H 
and AH). This work may be extended to solve the problem when length of the 
splitter plate will be more or less than the present case It may be possible to 
use sinuous trailing edge splitter plate Also the circular tube could be replaced 
by oval tube We can try to extend this to show the effects of splitter plate for 
very high Reynolds number We can aim for finding out the length of the splitter 
plate for which there would be maximum heat transfer corresponding to a given 
Reynolds number The effect of the length of a splitter plate on Strouhal number 
may be a interesting study to be investigated Here m the present computational 
study only a circular tube placed m a channel has been considered. This study 
can be extended for the tube banks 

A circular tube with an annular fin of outer diameter of 1 5D has been investi- 
gated. This work can be extended for the different outer diameter of the annular 
fin This study would be attractive if it is extended by attaching two or more fin 

m the prescribed geometry 
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